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ABSTRACT 


This paper presents two types of multi-item deterministic 
constrained inventory models which relate to the constraints 
of the equivalent of a direct support unit in the Federal 
German Army; namely, that there is only a limited amount of 
Space available for storage of items. One is a continuous 
review model with no stockouts allowed and the other is a 
periodic review model where backorders are allowed. In 
contrast to previous models analyzed by Naddor and others, 
time phasing of the arrivals of orders of different items is 
allowed. The benefit derived from this phasing is the 
reduction in the average amounts of on-hand inventory and 
backorders. Solutions are obtained for the two-item case 
and the n-item case is discussed. Examples illustrate the 


results. 
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I. INTRODUCTION 


Military forces have to be supplied with all kinds of 
items in order to complete their missions. Because many of 
the needed items cannot be purchased from outside suppliers 
at the time the forces need them, the items have to be 
purchased and stored earlier than actually neegea 

From the point of view of the ultimate user it would be 
best to store all such needed goods as near as possible to 
his own iocation, so that whenever an item is needed it can 
be delivered immediately. However, there are obviously 
constraints which make this course of action impossible. 

First, there is a budget constraint which allows only a 
certain amount of money to be spent on purchases of equipment, 
Spare parts and other needed goods and, second, there is a 
space constraint because the supply support units are required 
to be mobile so that they can stay in close proximity to the > 
combat units they support. In the case of direct support 
units for land forces the budget constraint is either non- 
existent or overwhelmed by the space constraint. 

In the Federal German Army these units are called BVP's 
and each consits of a fleet of trucks. They are equivalent 
to the Direct Support Units (DSU's) of the U. S. Army. The 
BVP's constitute the last echelon of the German Army's supply 
system and are not involved in buying goods from outside 


Suppliers. Therefore the budget constraint on purchases 1s 





nonexistent. However, the number of available trucks for 

a BVP is limited by its operational requirements and, to a 
certain extent, by the available budget for the armed 
services. Therefore any relevant inventory model has to take 
a space constraint into consideration when calculating optimal 
operating policies. 

Two models will be analyzed in this thesis. The first 
will be a multi-item deterministic lot size model subject to 
a space restriction which requires that the amount in 
inventory at any time cannot exceed the available storage 
Space and a restriction that the time between orders for all 
items must be the same. The variables which are being 
optimized in this system are the time between orders for each 
item and the times between orders for different items. The 
notion of allowing orders for different items to arrive at 
Peerent canes wlll be referred to as time phasing of the 
orders for different items. The solution for this model is 
obtained for the two-item case and the procedure for solving 
the n-item case is outlined. The second model is a multi-item 
deterministic order level system. It is subject to the same 
space and cycle time constraints as the lot size model. 
However, stockouts are allowed. The decision variables are 
the order levels of each of the items and the time phasing 
of the ordering of each item. This model is solved for the 


two-item case and discussed for the n-item case. 








mie PREVIOUS WORK IN THE AREA 


Pee CONSTRAINED LOT SIZE MODEL 

A multi-item lot size model with a single space constraint 
has been examined by Churchman [2] and Holt [6]. The model 
was defined in the usual form of deterministic inventory 
problems; that is, to minimize the total variable costs — 
reorder plus carrying costs — of the system. The space 
constraint required that the sum of the average space 
requirements for all items not exceed the total space availabie 
for these items. The assumptions which were made in the 
formulation of the problem were: 


1. ga Known constant demand rate; 
Z. infinite replenishing rate; 

3. constant lead time; 
4 


. shortages are not allowed. 

The cycle lengths of all items were not required to be 
identical. In addition, the time phasing of each item was 
not addressed. 

A disadvantage of this model is that only the sum of the 
average inventories must be no greater than the space available. 
This requirement can be satisfied mathematically and yet not 
be met in reality because the possible different cycle lengths 
may cause the order for more than one item to arrive at the 


same time, 





B. A CONSTRAINED ORDER LEVEL SYSTEM 

A multi-item deterministic order level model has been 
examined by Naddor [7]. This model used the first three 
assumptions of the above model by Holt and, in addition, 
required a fixed review interval for all items. Backorders 
were allowed however. The objective function for minimization 
was the total variable cost — carrying plus backorder costs — 
and the constraint required that the summation of the space 
requirements of the optimum order levels not exceed the total 
Space available. 

The major difference between Naddor's model and the one 
to be presented in this paper is that in his model all orders 


“js 


arrive at the same time, while the model being developed here 


a 


Lt 


Ww. 


‘ermits time phasing of arrivals of orders of different items. 


3 
His solution for the optimal order level for the n-item case 


is: 
a ~ Fy - a | 
ial Zal 

where 

Says order level of item i when lead time is zero; 

q; = replenishment size of item i; 

cj, = the carrying cost of item i per unit of time; 

Co; = the shortage cost of item i per unit of time; 

b; = unit space required for item 1; 

B = total warehouse space available; 
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and 


Ve. 





LE 


Peecelote en eOeMULTI-ITEM LOT SIZE SYSTEM 


The model which is developed in this chapter is ann item 


lot size model constrained by the amount of space which is 


available 


are made: 


for storage of inventory. The following assumptions 


[veCTanas bor Litem 1 LS Geterministic at a constant 
positive rate of a units per unit time. 


2. No shortages are allowed and the stock of an item will 
be immediately replenished whenever its inventory level 
reaches zero. 


Be. The 
4. The 
5. One 


6. The 


lead time 1S constant. 
replenishment rate is infinite. 
type of storage is usable by all items. 


Mivemeoryecarrvying COoSt 1S proportional to the 


average amount in inventory at any time; its dimension 
is $ per unit-time. 


7. ###The 


order cost is a constant which 1s incurred each 


time an order is placed; its dimension is $. 


The variables which will be optimized to produce minimum 


cost are the length T of the reorder period for every item 


and the time x at which an order for the ith item arrives 


after arrival of an order of item number 1 (x, = 0). The cost 


equation, 


as a function of the scheduling period for any 


Single item under these assumptions, is given by 


= wae 2 UE il 
C; (T) C,3r, 7/2 + c4,/1 (1) 
where Ci; = the inventory carrying cost, and 

C3, = the unit reorder cost. 


12 





This notation is that of Ref. 7, page 58. For n items, the 


total cost can be expressed as 


(Cc, 57; 


{7/2 + C3,/T) P 


given that all items are using the same reorder period length. 


The space constraint can be formulated as 


n 
ee Oren (tC). -<  B hOmcas 1 tt > FO, (2) 
al ii ~~ a 
1=1 
where 
B = total space available; 
b. = the space occupied by a unit of the ith item; 


T, (t)= the inventory level of the ith item at time t. 


Since maximum inventory levels occur at the times at which 
orders are received, these times correspond to critical points 
relative to the constraint. Therefore the constraint > 
only to be examined at these arrival times. This will result 
in n constraints for an n-item case because of the possibility 
that each item may arrive at a different time. 


To solve the problem let hs represent the ith Lagrange 


multiplier. Then the Lagrangian becomes 


n n 
Boas Eee. (T) + y dF, (3) 


ns 





where C; (T) is given by equation (1) and FE. represents the 
fem constraint function (describing (2) at each Xs 
i=1,2,...,n}). The solution to equation (3) will be presented 
for cases which have two and n items. The two-item case is 


presented to illustrate the general solution procedure. 
Pee OO LUTION OF THE TWO-ITEM CASE 
Equation (3) for the two-item case reduces to: 


L = C,45\7/2 ~ Ca,/T +C)5r5T/2 ~ C35/T + VF, + AUWF 


a 

Typical fluctuations of the inventory levels of the items can 
be represented as shown in Figure l. 

From Figure 2 it can be seen that the constraint function 


at time zero when item number 1 arrives will be 


¢ 


ana thatthe constraint function at time Xo when item number 


2 arrives will be 


byr,(T ~ Xo) + boroT oe: 


Conversion of inequalities into equalities for use in the 


Lagrangian yields: 


2 : 
qi b,r,T + Dor5X5 BS V4 - B , (4) 


ry 
It 


2 p 
gXoT + V5 - B ji (5) 


| 


Fe = by r,(T - Xo) + b 


where Vy and Vo are slack variables. 





item no. | 


mem no. 2 





FIGURE 1: Fluctuation of inventory levels for two items for 
the deterministic lot size model. 


LS 





-—4 
-—— 

ct 
al 


occupted space 


byryT = boro! 


DaPyX_ = by (T-x>) 





FIGURE 2: Space occupied by two items when bor, = bol, 
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Comormang (3), (4) and (5) gives: 


L = C44%,T/2 + C3,/T 4 Cyor5T/2 + C35/T 
Tevet Yeuee+ DbAr.x. + V a B) 
ee = 232-2 1 
2 
+ Ao (by xr, (T Xo) + boroT + V5 - B) . (6) 


Taking partial derivatives of the Lagrangian with respect to 
More Vo Xor AR dy and hos setting them equal to zero, and 


Simplifying, yields: 


AV) = 0 ; (7) 
hoVo = Q : (8) 
dybor, hob ry ae) ; (9) 


Ye 2 
C1423 /2 - C3,/T + Cy9%5/2 C35/T + hybi ry 


+ dobry + Ao bor, = 0 ; (10) 
2 S : 
b,r,T + bjr5X. + V1 =—~ea = 6 ; (11) 
2 2 
b,r,(T- x) + b5roT + V5 - B = 0 : (2) 


Equations (7), (8) and (9) indicate that there are only 


two possible cases: 


17 





¥ 0 and ho *# O which implies that Vi = V5 = 0; 


a. AS = ho = 0 , and the V; values can be arbitrary. 
When hi * 0 and ho # 0, both constraints are binding. Solving 
equations (ll) and (12) with Vy = V5 = 0 for T and Xo 
respectively, yields: 

2 Ih ai 


l +b rj bor, + bo ro ys (13) 


il 


2 
* 
T B(b xr, +b5r,)/(b, “x 


a 2: eee 202 
Bboro/ (by ry + b,x boro + bo r5 ao (14) 


Substituting T and Xo BiEomtheuretal cost equation yields 


the minimum cost Cy where 
C = B(c 4%) +f ,9%5) (bir) + bor) /2D + D(c3, +C35)/A ; 


and ; 


me. 2 2 
D= by ry + bir bor, =F b. ro : 
A B(bir, + bor.) 
The ratio 
* 

2G) Dory 

—- = 

Ab | Dir, + bor. 


18 





provides a convenient means of comparing x5, to T for this 


completely constrained case. If bir, - bors then 


* * 
ae Ee chen 2 yl te 15 ys 
a a Lh eS: e = > ° 171 5X5 


For the unconstrained case, equation (9) requires that 
if one of the A values equals zero, the other has to be zero 
too. The optimal unconstrained value of T, denoted as ye 


can be determined from equation (10). It is 


(15) 





Equation (15) is the same expression as derived by Naddor [7]. 


The value for x. 1s dependent, however, on the value of B. 


Z 


The value of B which corresponds to the dividing point between 
cases 1 and 2 is obtained by equating equations (13) and (15) 
and solving for B. This value is denoted as By and is given 


by equation (16) 
a D 


When B SB) , equations (13) and (14) give the values of 7 


* x 
and x, . When B > B = Th and the value of Xo is no 


2 lias? 
longer restricted to that given by (14). 


19 





If B 2 Bo where Bo = (bj xr, +bor,) ste pee eEnen Ko can be 
any value in the interval 0 SX < T,- This corresponds to 
B being so large that even if orders for both items arrive 
Simultaneously they will not be faced with a storage 
restriction. When By <0 8) << Bo , the storage constraint will 
not be restrictive provided that Xo is so selected that the 
constraint is not exceeded at the time either item arrives. 

When the first item arrives (at t = 0) the total combined 
on-hand inventory consumes birt, a boroX5 amount of 
storage space. When the second item arrives (at t = Xo) the 
total combined on-hand inventory consumes space in the amount 


of bir, (T, - X5) + b The value of Xo must be selected 


uaa ‘ 
so that the total combined on-hand inventory at either time 


does not exceed B. Thus if Xo lies in the interval 








, (17) 


then both conditions will be satisfied when By < B< B,. 
When B = Bo the interval given by (17) simplifies to 
0 S X5 st - As B decreases below Bor the upper bound 
decreases and the lower bound increases at the same rate. 
The interval reduces to a single value for Xx, at B= Bi; 
the corresponding Xo value is given by equation (14). Asa 
consequence of this behavior of the interval for B between 
By and Bor it follows that equation (14) gives an optimal 
fom@eall B > B 


value of x , and the only optimal value for 


2 a 
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Bee SOLUETON OF THE n-ITEM CASE 


The Lagrangian for n items is: 


1g) 
C,3,r,T/2 + C,,/T) aus ne A.FP, / (18) 


with the constraint functions given below. These functions 


assume 0 < X5 < X23 S$ se S$ xX, . 


Ze 
Ze 

r, (T-x,) + bor r5T +b 3r 3 (X37 X) + ok +b r ne y7*) + V5 = 8B 
= — 2 _ 

bj x, (T-x,) + box 7 (T- (x5 x5) ) “ie bor.T + oon + Pann x3) a = 8B 
Ze 

b 171 (TX, a bor 5 (T - a y7%)) + oe + bor (47x) enc 5B 
| 5 

bjr, (I+) + bor, (T - (x, -%) + oe + borat + ve =a 


After forming the Lagrangian, taking partial derivatives with 


respect to Xo1 Xgr veer X and Vye Vos eee and setting 


n 6 a 
them equal to zero, it can be shown that either all of the 
Lagrange multipliers are zero or all slack variables are zero. 


The partial derivatives with respect to Aye hos ater hve 


2h 





when the V,'S are all zero give back equality constraints at 


t = 0X5 nee eX, - Solving these equations for T, Xor seen Xe 


n 
yields: 
n 
Bee bor 
ko j=] 12 
T -_ 7p) ee f (19) 
a 
Be >) - Mb. r 
. ees 
7 = > © r La 1 eg Ae | 7 (20) 
where 
n > 2 Deed n 
D = 16 aa: a Y bx, » Diese. . (221) 
= i=l j=itl JJ 


The resulting value of Cor the minimum total variable cost, 


Tht 


n 
a 


c2./B Weer. +...tber.) ; 
eT ea 31 re Taehle 


n 
Co = B {2 C)505/2 (b)r,+...tb xr) /D + D ( 


In the unconstrained case, optimal T is obtained from 





(22) 
DT, 
B) = = ae (23) 
a oe 


lag 








n 
Be = temo. ©. T ; (24) 


# 
And, if B > Boy Enenes, “= Pe and each x; Can take on any 


Value in the interval 0 < x. < sae - When B 


eS Bau< Bo then 
the x, 's must be chosen in away that the constraints are 
not violated. As B decreases below Bos the upper bound 
decreases and the lower bound increases. However, equation 
(20) has been obtained for the case where the constraints 
have to be met, 1.e. where we had equalities in all of the 
constraint functions. Changing the equalities to inequalities 
will result in an infinite set of optimal feasible solutions. 
But x; as given by equation (20) will always be in this set 
of optimal feasible solutions because the solution to the 
more restricted problem will always be an answer to the less 
restricted problem. 
The set of optimal solutions creates the ranges for the 
xX values. They could be obtained using linear programming 
in the following way: 
1. Possible objective functions are 
a) min Xe imM—9 255,252,708 to get the smallest lower 
bound for x7 
b) max Xa (ie——e oye a he cOrgeta cae Greatest upper 
moume L£Ox Xs 


2. The constraints are all of those shown below 


equation (18). 
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For solving problems which have n > 2 items the flow chart 


of Figure 3 is helpful and summarizes the preceding analysis. 


C. EXAMPLES 


Two examples are presented below to illustrate the model. 


Example 1 


Consider a warehouse which stores two different items 


whose parameters are: 


meem 1: Item 2: 

mr, = 200 units per year ro = 250 units per year 
b, = 5 £t°/unit b, = 3 £t°/unit 

C11 = $ 2 per unit-year Cj5= $ l per unit-year 
C3, = S30 Cz5= 9 25 


If total available warehouse space is 400 cubic feet, what 
* ‘ 
Sa-mencwOpeimal Order period fT for each item and the optimal 
* 
arrival time Xo for item number 2 after an order of item 


number 1 arrives? 


Calculating Th and B, yields: 


a 


2 ; _ 3 
a = 0.41 years ; aS 543 £t : 


Because B < B, , both storage constraints are binding. 
* * 
Therefore we get T = 0.30 years and Xo = 0.13 years by 


using equations (13) and (14). 
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Calculate: 
T= equation (22) 
B, = equation 28) 


T* = equation (19} 
J ox 
x. = equation (20) 


calculate: 
equation (24) 





yes 





x e 
Ae equation (20) 





FIGURE 3: Flowchart for solving lot size model problems 
having n > 2 items. 


iD) 





Example 2 
Suppose the data of example number 1 remains unchanged 


except for B. Assume B is now 600 ft, 


The values for Ty and Bi are unchanged. But since B is now 
greater than By. evaluation of Bo is required. It is 
calculated to be 719 ft, 

Because By <B < Bo ; T = Th = 0.41 years and Xo lies 


in the interval given by (17); that is, 
Diz years ees Xo < 0.25 years . 


Checking the feasibility at time t = 0 shows that, using the 


* 
lower bound as the x, value, the occupied warehouse space 


2 
* 
would be 500 eee using the upper bound as the Xo value, 
600 ft would be needed. The respective space needs for 
* 
time t = x. would be 600 ft” and 467 ft>. 


Z ‘ 


D. COMMENTS 
The constrained multi-item lot size system developed in 

this chapter differs from that of Holt [6] in the following 

ways: 

1. Holt's model allows a different cycle length for each 
item, while a common cycle length for all products is 
used in the above model. 

2. Holt's constraint required that the sum of average 


inventories be no greater than the space available, while 


26 





the constraint function used in the above model requires 
that the space occupied by the inventory at any time be 
no greater than the total space available. 
3. The time phasing between orders for different items is 
calculated in the above model, while the time at which 
an order arrives is not regulated by Holt. 
The second aspect 1S perhaps the most important difference 
because the actual space restriction in Holt's model is a 
nebulous notion. Realistically, the space constraint restricts 
the maximum on-hand inventory rather than the average level. 
It is interesting to note that the solution for the 
scheduling period can also be used to calculate a minimum 
amount of storage space needed given a specific scheduling 
period. The equation for the warehouse space in terms of the 


scheduling period is 


where D is given by equation (21). 
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LV. 


PVCONSTRAINED MULTI-ITEM ORDER LEVEL SYSTEM 


A. THE MODEL 


The system which is analyzed in this chapter is an n-item 


order 


level or periodic review system which is constrained by 


the total amount of space which 1s available for storage of 


the items. The analysis of this system will be conducted 


under 


1 


2. 


2) 


the following assumptions: 


Demand for item 11S deterministic at a constant rate 
of rs units per time. 


The period between replenishments or scheduling period 
is a fixed constant T. 


The replenishment quantity of the ith item raises the 
inventory at the time the quantity is received to the 
order level S; after backorders are made up. 

The replenishment rate 1s infinite. 

The lead time is constant. 

The same type of storage is used by all items. 

The inventory carrying cost is proportional to the ‘ 
average amount in inventory. Its dimension is $ per 


mnit time. 


The backorder cost is proportional to the average shortage 
level. Its dimension is $ per unit time. 


The order costs are ignored because they are a fixed cost 
iasche scheduling period is a prescribed constant. 


From these assumptions it can be seen that the lot sizes which 


are being considered are constant since 


where 


25 a as in 


q; represents the ith item lot size. 
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The variables which are to be optimized in this system 
are the order levels (S; for the ith item) and the points in 


time at which each item arrives (x, for the ith item, where 


The cost for one scheduling period for the ith item as a 


function of the order level can be shown to be [7]: 


2 
C1383 Co; (43 7S8;) 


2 Er Plow 
1 1 


inventory carrying cost, 


= 
= 
49) 
= 
(D 
Q 
— 
KF 
il 


inventory shortage cost (using Naddor's 


and Co; 
notation). 


Then the total cost equation for an n-item model can be 


expressed as: 


2 2 
nN Cc Sir Cc . (7 e=S.) 
= Jal ak Zila gee 8 
C(S)1So1-++18)) 7 2 (37-7 «Ft Ze. ) 
i=l a ale 


The storage constraint can be formulated as: 


n 

Tee, (ct) <8 fOorwalt t€ = 0 , 

: ee a ~ 

i=l 
where B = total space available, 

b. = amount of space which is occupied by one unit 
a Of cele, 
I; (t) = inventory level of the ith item at any point 


mae time t. 
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As in the preceding model, the maximum amount in inventory 

of any item occurs immediately after the ordered item arrives. 
m@enrerore the arrival time is a critical time with regard to 
the space constraint. Because a constraint will be required 
for the time of arrival of each item, there will be n 
constraints for n items. 

Figure 4 illustrates the two possibilities in fluctuations 
of inventory for two items. Two cases are possible because, 
at the end of each scheduling period, an order for one item 
arrives and the other item is either on backorder or in stock. 
Care must be taken in formulating the constraints because when 
an item is on backorder its inventory level is negative but 


the amount of its occupied space is zero. 


Seo OLULIONS FOR THE TWO-ITEM CASE 
The problem for the two-item case can be stated as 


follows: 


Minimize 


2 a 
- ei | Co, TyT Sy) Aaa oan ae 
rer. mr new! = “7. ra ahs P) 
2rjt rT 2roi aro 
subject to the constraints 
bjS) + max[0, bo (S, ~ xr, (T-x,))] ae 5} : (26) 
b5S5 + max([0, b, (8) - 1 %5) I <a : (27) 
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Case | 





Case 2 





FIGURE 4: Fluctuation of inventory levels for two items for 
the deterministic order level model. 





Constraint (26) reflects the combination of the constraint 
cases at t = 0. Whenever Sy — ro (T-x,) < 0 , which means that 
item number 2 is on backorder, this constraint reduces to the 
mecstriction that bjS, be less than or equal to B. Similarly, 
constraint (27) reflects the combination of cases at t = Xo 

ie totally Unconstrained Case 

Tf the problem is unconstrained, then the optimal value 
EOr Sy and S5 can be obtained by simply taking the partial 
derivatives of equation (25) with respect to Sy and S., setting 
them equal to zero and solving for Sy ands co. The solutions; 


2 


denoted as S and S 7, are: 
EV u2 


Be = 2a (28) 
ul C11 $251 
Cee 
one) a (29) 
ee 22 


which are the same expressions as in Naddor's unconstrained 
> + = 

ease [7/7]. Because B > biSi1 b5S..5 By , even simultaneous 

arrivals are not Space constrained and x, Can assume any 

meme in the interval O < Xo <T. 


2: Bete ele, Unconstrained Case 


When B < B then simultaneous arrivals would 


3 a 
violate the space constraint. If B is not too restrictive, 
the problem could still be unrestricted relative to Sul and/or 


S42 if careful time phasing of the orders is done. 
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The individual constraints were: 


b,S) + max[0, bo (So — rg (T-x,))] 


A 
w 


boSo 


A 
e0) 


+ max[0, b,(S,;- 4 X5) ] 


Notice that Sul and S32 will not necessarily violate the 


Space constraint if B > Bo and Xo 1s chosen carefully, where 


Bo = max{b,S.4. boS..5] a (30) 


They will therefore remain optimal. 


If we assume that both pb Sy < Band bos < B 


Leauge u2 


then the functions max(0O, ba(S.5 7 Fg (T-x,))] and 

max[0, bi (S44 7 1%) ] at t = 0 and t = X51 respectively, 
could be positive. At t = 0 the value of the max function 
1s not to exceed the total space available minus the space 


occupied by item number 1. A feasible positive max function 


at t = 0 satisfies: 





0 < bo (S45 — Fo (T-x,)) [R= bis : 
@e, in terms of Xo 
Ss S B-b.S 
ae We - tee oar 
2 2 252 


Sie) 





S 
u2 : 
5 Saal ca then bo (So r(T-x.) ) 0 eee I Xo is 


greater than the right-hand bound then the total inventory 


ie xX 


at t = 0 would exceed B, the storage restriction value. 

When t = Xo the value of the positive max function 
is not to exceed the total space available minus the space 
occupied by the optimum order level of item number 2. This 
means that 


Cem onie a>) = BP 5S5o 


or, in terms of Xo 


nn, Gl aan < : (32) 





S 
; ae 
Bietnis Case 1f 8 Te , then b, (S474 1%) S 0 


and if Xo < left-hand bound, the total on-hand inventory 


at t = x, would exceed the storage restriction. 


2 


In comparing inequalities (31) and (32) it is clear 


that the space constraint will be only met if: 


Sam © -o°yo Seg 1 aa 
wee. tC Bee 
1 at, => Dae, 


However, feasibility requires that the left-hand side of (33) 


be less than or equal to the right-hand side. 
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The two bounds of (33) are equal when B = By where 


Bees °5 > Bo (34) 
and 
By = 548.7 * 228.2 | (35) 
Os ead on ay 
Bo = [5 : ct - JT . (36) 
17, toro 


Therefore, if Bo < By. GB =< B. then the problem is 

* 
unconstrained with respect to Sy and So provided Xo is chosen 
within the boundaries of inequality (33). 


Simplifying (33) results in: 


B,-B B,-B8 


It follows that when B = B the range given by (37) is 


3° 
O< Xo < Ti when B < B. the range becomes narrower with the 
upper bound decreasing at a rate 1/boro5 and the lower bound 
increasing at a rate -l/b,r, per one unit decrease in the 


value of B. We note in passing that if B= bis = boS 19 j 


the range (33) reduces to 
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es ; : ; 
mye condition resulting in Be max[b,S.4,b5S_ 5] 1S a case 
‘ where both constraints are binding for any B < By while the 
max functions are both still positive for B near Bye 

3. Constrained Cases 


aA. Bos and B 


x 2 0 
RIGE Bes > max[b,S 1 +b9S\.5] and B< By. the 


< B< B : 
5S 


constraints when the max functions are both positive are: 


[A 
ee) 


1A 
eo) 


Subtracting the second inequality from the first gives: 





Da B 
~ Zee 0 
x = | | es : (38) 
2 byr, t bor, bir, 
* 
Thus, Xo 1s independent of B provided that both max functions 
* ¢ 
remain positive. Substitution of Xo into either constraint 
gives: 
54S) + b5So ee + br xX F 


[A 
es) 
+ 
to 


Therefore we have the very simple relationship that, whatever 


values of Si or Sj, are chosen, we must have 


2 
; * * 
b,S, + b5S. 6 Bo : (39) 


36 











Now inequality (39) does not reflect the fact that the space 
constraints at t = 0 and t= Xo would be violated if either 

* * 
bjS, a Or boS, mebawee Coererore, the following conditions 


must also be imposed: 


[A 
wo 


(40) 


a) 


oy 
Ww 
1A 
es, 


The form of the resulting constraint region is illustrated 
in Figure 5(a). The unconstrained optimal solution (SS y2) 
lies outside this region. Because the objective function is 
convex in Sy and S. and the constraint region is convex, it 
follows that the constrained optimum lies on the boundary of 
the feasible region. 


If B is reduced to B, we get the condition illustrated 


0 
in Figure 5(b) and when B < By we get the result of Figure 
5(c). These cases will be analyzed in the next section. 

To search for the optimal solution when Bo < B < Bes 
we evaluate Sy and S5 using only one segment of the boundary 
at a time. We begin by assuming that the constrained Si and 
S. lie on the line given by (39) with equality. We can then 
use (39) to solve for S5 as a function of Sy: substitute that 
function into the objective function, differentiate with 


respect to Sy: set the derivative to zero and solve for Sy: 


We get 


37) 





FIGURES: 


ie 


© = unconstrained solution 






Feasible regions for the two-item deterministic 
order level mode) where aro 7 paoe Ss By ; 
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0 1 22 
oe pe) tI) See 
S, = 5, = ; (12) 
2a 1 ey: 
[=] + [pl ie 
ul u2 
aK B+B, b) 4 
S5 = So =a | By = Sy (43) 


If either S, Sg S. violates the DiS) Seki eke boS5 < B 
constraints then the optimal constrained values lie on either 
of these latter constraint boundaries. 

iets S, is larger than B/by then the largest feasible 
Sy value is B/b, - This then will be the optimal value of Si: 


* 
Substitution of oae = B/b, into (igmobjective function and 


v 


Zt 


differentiating with respect to S5 will give 32 when the 
PS eedulvewis Se6t eto zero. “iherefore, optimal constrained _ 


S5 is the largest value the constraints will allow without 
* 


mecucCIng Si: This corresponds to the S5 value at the inter- 
* 
Section of (39) and (40); that is S5 = By/b.- 
* 
Similar arguments lead to Sy = By/by and 
* A 
So = B/b. when So is greater than B/b. 
* 
In summary, 1f By < B< By then Xo is given by 
* * 
equation (38). To determine Sy and So: 
Pe Calculate Sy: S5 


eet eo eey bays, < B/b, then S, =S,. 8S, = 8, 


oe 





Dee Bis Bo and B < By : 
Figure 5(b) and 5(c) emphasize that when B < Bo 
the constraint given by (35) lies outside the boundaries of 
the feasible region and is trivially satisfied by any solution 
in the feasible region. It follows from the arguments in the 


* * 
preceding section (a) that Sy = B/by and S. = B/b. and 


2 
we realize immediately that both max functions must be zero. 
The analysis of the preceding section showed that 


if 


then both max functions will remain zero. Substitution of 


* * 
Si and S5 into this range gives 


ee 








This range reduces to = By/D4ry (equation (38)) when 


ay: 
B= Bo: the range increases to 0 < X5 < T as B is reduced 
Io AS Oe 

Figure 6 summarizes the ranges of optimal Xn des 


the case where B, > Bo . From this figure it becomes obvious 


that x, = B/D ry is optimal regardless of the value of B. 
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FIGURE 6: The range of optimal Xo aS a function of B 
for the case where BY > Bo : 


4 


a 





cer B.A 7 B and B< B 
a 2 


In the constrained case where Bo > By. pur ener 


reduction into cases where DiS zs DoS 5 P DoS 15 > D541 


and bu = 545.49 are necessary. We will also make use of 


By = min{b,S,4/b55,,5] : CaS) 


(ly): Be = biSu. By = b and By < B<B 


Be)! 2 


When Baws B < B tien te £LirSt constraint 


mompecome binding 1s that at t = 0. Its max function will be 
zero. An appreciation of what 1s taking place can be obtained 


by looking at the x, range given by (33). When Bote Be then 


2 


the upper and lower bounds do not reduce to the same value 


(at B = BL) before the B-b Sul term in the upper bound 


i 


goes to zero. The resulting range for X5 from (33) is now: 





S 
<T- —— (46) ' 


* 
when Si = B/by 1s substituted for Sut The lower bound 1s 
Suly valid however as long as B > By = boS.9 wee The Opermal 


value of S. is § 


Z u2~- 


Inspection of (46) shows that it is 





independent of B provided that B, <B<B, . When the value 
of B is reduced to Bue (46) reduces to (44); that is 
B 
<< < T- —_—_— _ , 
biry —-  2—- bir: 
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ena, since S. must be B/b. FOnevalil B \< By , this range 
must be valid for all such B. 

The analysis in section b above showed that 
(44) gave feasible Xo values sas tong as B < Bo - Tne analysis 
in this section then suggests that B, must be less than or 
equal to Bo 1£f (46) and (44) are to be feasible and identical 


at B= B Now the range of (46) is feasible if 


1° 





Collecting like terms results in 


lola 

ilies 1: 

bes < { -—_————- ] T ' 
Dei br, tbor, 


which reduces to By < Bo 
In summary, 1f B. = bjiS4 = 8B and 


* 
B, = boS.5 <B < Bg then x, is any value in the range given 


* e 
by (46) and oa B/b, 1 So = S42 . When B < B, we again 


* 
have x 


* * 
> given by (44) and = = B/b, 1» Sy = B/b, 


Figure 7 shows the range of optimal X5 as 
Benicio or 8B £or this case. It is immediately evident 
that x, = Bo/byry is optimal regardless of the value of B. 
and B 


(Zi b= =" b.S , B, = Db 


2 ene 1 L?ul it 


By similar arguments to those given above in 


< 2B V< Bo 


* * 
section (1) we get S, = Sy and S5 = B/b. - The x, range 


from (33) reduces to 
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FIGURE 7: The range of optimal Xo aS a function of B 
for the case where BY < B, and By < Bo . 
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oo < x, Se ie — ’ (47) 
IL QT 
* 
when By <b. Bo 5 illgitets\ Jes By ; Sy = B/b, ; 
* * 
S5 = B/b. and xX» can be obtained from (44). Figure 7 also 


illustrates the range of optimal x», values for this case. 


1 ~ PySu1 = P2oSy2 = Bo 


In this case, the interval By < B<B 


(oy <B 


* 
does not exist and therefore we have Sy = Sul P S5 = Su2 
and the Xo interval is given by (37) when van o>) 2 Ba . 


* * 
then §& SS... =UB/b and the x 


1 1 i Z 2 2 


interval is given by (44). Figure 8 shows the optimal range 


Wen B < B = B/b 
of X5 as a Piet kOon Of Bb. 
Having now examined all possible cases it can be 


seen that there are three major possibilities: 


Meese. the problem 1s unconstrained if B > By and 
partially unconstrained (the range of Xo is limited) 
* * 
3° SOlMmETons. 1oOr Sy and S5 


are the economic order levels S and S .. 
ul u2 


if max [B. By] “Sie Se! 


Second: If the storage constraint is between min[B, Bo] and 
max{B. ,B,] then the solutions depend on the relation- 
ship of the parameters to each other and to B. They 
can be classified according to two ranges for B: 


namely, Bo < B< By. and By <B< Bo - 
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FIGURE 8: The range of optimal Xo as a fUnGHlon on 6 
for the case where ER < Bo and EB, = Bo 
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a. * 
Solutions for Sy and So are: 


* a * 
ie Sy aed Sy c So aa 
* * 
Be Sy = B/b, ; So 
* * 
* * 
* * 
Dis Sy = B/by ; So = 


The first three solutions correspond 


* 
case and X5 can only have one value. 


correspond to the B, <B < B, case 


IL 


values within a computed interval. 


Taxed: In all cases where B Suusiol 


Constraint 21S SO reStrictive 
be stored at a time. Hence 


ane the interval for Xo 


GOmrEnNe 1B. < B < by. 


0 


The last two solutions 


* 
and X5 Can take on any 


[B,1Bo ] the space 


that only one item can 


* * 
S, = B/by and So = B/Do51 


J 


expands as B decreases. 


The details of the optimal solutions of the two-item case are 


summarized by the flowchart presented in Figure 9. 


eee Oa LTEM CASE 


The n-item case can be easily formulated using the same 


approach as for the two-item case. 


The general n constraint 


equations can be written for the’times at which the items 


arrive uSing the max functions. 
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Calculate: 


= equation (28) 


= equation (29) 


equation (35) 


Calculate: 


equation (36) 


30) 


( 
equation (34) 
equation ( 

( 


equation (45) 










* 
yes | = uy) 
2 = 
B > max{B,,B,) § 35 = Sy9 
* 
xX, in range 33) 


no 


* 
FIGURE 9: Flowchart for determining ae - and x, for the 
two-item deterministic order level model. 
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S 
yes 
s 


X 


* 
] 
= os 
a 
* 
> in range (47) 





no 





FIGURE 9: (continuation) 
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Calculate: 


equation (42) 


equation (43) 





FIGURE 9: (continuation) 
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ire principle difficulty involved in trying to solve this 
problem is the number of possible combinations which must be 
considered because of the various max functions. Because of 
its complexities the n-item case will not be addressed in 


this thesis. 


D. EXAMPLES 


es Example Number l1 


A warehouse stores two items with the following 


paemameters : 
Item Number l1: PEcnhmNumoer 2° 
r, = 200 units per year; i, = Z50MUNLES Per year; 
eee Zune, by = |S Becta at 
C11 = $ 2 per unit-year; Cio = S$ 1 per unit-year; 
Co, = $30 per unit-year; Coo = $25 per unit-year. 


Total warehouse space available is 600 oe the scheduling 


period T is one year. How many units of each item should be 


ordered and when should the orders arrive? 


Following the steps of Figure 9, we first calculate 


Sul! Su2 and B3. 


egies 3 ~S = 240.4 ; B,= 1658 . 
u 3 


ou 





Since B= 600 < B., , we must calculate Bor De By and Bo. 
BA = 429° B. =m s § Bs = J2i -; Bo = 937 


We see that By. > Bo ana. B. < B and we will have to 


calculate Sy and So - 


FOS. 6 


WY > 
I 
ps 
© 
Le) 
has 
WK 
it 


Now, because S, < ~=120 and S, < = 200, the 

a by Zz bo 
optimal values for Si and Sj are Sy and S5- The corresponding 
Xo value 1s given by equation (38). In summary, the optimal 


solutions are: 


Cp) 
il 


* 
106.4 units ; So orion ts, 


0.429 years. 


Total costs TC, uSing equation (25), are $ 1125 per year. 


2. Example Number 2 


The available warehouse space is the same as in 


Example number 1. However, the parameters of the items are: 


a2 





Item Number l1: Item Number 2: 


es 200 units per year; ro = 250 units per year; 
b, = 5 £t>/unit; by = 3 £t°/unit; 

C11 = $ 5 per unit-year; Cio = S$ l per unit-year; 
Co) = $100 per unit-year; Con = $ 3 per unit-year. 


=—~1907; Ss Sono eb. — 1512 


Since B = 600 is less than B. calculations of Bor Bye Bir 


3 


and Bo are required. 


By at 29. Dae OS 5B. —-562.5 4 By = 950 


Because B_ > Bo and Bo < B we have to calculate Sy and 


So @ 


> 


Because §5 > BL 
l Dy 


The optimal feasible solutions are: 


Il 


* * 
Sy B/D, = 120 units ; S, = By/bo = 143 units ; 


+ 


= By/b ry = 0.429 years. 


LOetimCosts musing equation (25), are $ 1890 per year. 


ope 


= 120 , Sy 1s not a feasible value for S3- 





E. COMMENTS 

The major difference between the model examined by 
Naddor and the model just developed is that the different 
items in Naddor's model all arrive at the same time, while 
the model being considered in this paper permits time phasing 
of the arrival time for each item. This time phasing allows 
a higher average inventory level for each item and, conse- 
quently, a reduced average shortage level and a lower total 
average annual cost for the system. For the data of Example 
number one the optimal order levels using Naddor's model are 
Sy = 51.3 units per year and an 114.5 units per year, 
yielding total average annual costs of $2600 (using Equation 
(25)). Time phasing therefore "saves" $1475 per year. 

In addition to Naddor's model being more expensive, the 
solutions are not always feasible. Naddor's analysis 
implies that a negative optimal order level "creates" 
additional available space. In the model of this paper this 
problem has been carefully considered in the max functions. 
Naddor does not address this problem at all: Using Naddor's 
results to solve Example number 2 we get S, = i797 Unlts and 
So = -94 units. Clearly, DiS) > B and hence Si = 177 units 


EomePatecasiie le. 
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V. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

Chapter III provides a minimum cost (ordering plus 
holding) solution to a multi-item lot size inventory problem 
Subject to a constraint on the maximum total space available 
for storing the inventory. The variables which were opti- 
mized were the scheduling period (time between orders of a 
Single item) and the inter-order times between the different 
items. The assumptions which were made in the formulation 
of this model were a deterministic demand, infinite replen- 
ishment rate, no shortages allowed and a common scheduling 
period for all items. 

Chapter IV provides a minimum cost (holding plus back- 
order) solution for a two-item system under periodic review 
and subject to the same space constraint as in Chapter III. 
The variables which were optimized were the order levels 
for the two items and the inter-order time between them. 

The solutions have several possible forms depending upon 

the interactions of the constraint and the parameters of 

the items. The assumptions which were made in the formula- 
tion of the problem were a deterministic demand, infinite 
replenishment rate and a prescribed interval between reviews 


which is identical for all items. 
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B. RECOMMENDATIONS 

Basic to the models presented in this thesis is the 
assumption of deterministic demand. In reality, we are 
faced with random demand patterns and therefore it is impossi- 
ble to determine with certainty the future states of an inven- 
tory system. However, it has been shown that examination of 
deterministic inventory models can provide a better under- 
standing of the situation being modelled and may suggest 
approximations which would apply in the case of random 
demand [4]. This thesis provides such an understanding and 
one of the next steps towards improving existing severely 
constrained multi-item systems is to obtain demand data and 
classifications of items as either needing continuous 
reviewing or needing only a periodic review. 

The analysis of the lot size model in Chapter III was 
conducted because there are a certain number of items which 
Spommcenever be out of stock in a BVP. Further research 
should be able to provide a listing of these items anda 
probability distribution to describe their underlying demand. 
Because these items should be continuously reviewed to avoid 
stockouts they should probably be stored together on one or 
several trucks to facilitate careful monitoring. 

The periodic review model of Chapter IV was presented 
because the majority of the items which have to be stored by 
a BVP are not of such vital importance that they cannot be 


on backorder a certain amount of time. This model emphasized 
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the importance of knowing the values of time-dependent 
backorder costs. Such costs are needed and military essen- 
tiality coding should be a valuable help in appraising such 
costs for each item, thus insuring that the more important 
items have higher backorder costs and will therefore be less 
on backorder in an economic inventory system. The demand 
distribution will also be needed before the model results 
can be extended to a random demand situation. 

Though the problems of solving the n-item case in 
Chapter IV seem to be severe for large n, situations involving 
large n valves probably do not occur ina BVP. As has been 
already mentioned, a BVP consists of a fleet of trucks. 
If the items which are to be stored on one truck are selected 
in a way that all items have similar parameters — low cost 
and high demand for example — then it could be possible to 
form classes of items and treat each class as one item. Thus 
the classes per truck could perhaps be reduced to two or 
three. Grouping of items into classes of similar demand 
would also eliminate the need to determine a probability 
distribution for each item; the demand distribution for the 
class would suffice. Of course the notion of "similar demand" 
is rather nebulous and would require some further defining 
before grouping could be made. 

While demand and cost data are being gathered, further 
analytical work should proceed using the ideas behind the 


approximate random demand models of Chapter 4 and 5 of Hadley 
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eiomwoitan [4]. En particular the (Q, r) model of Chapter 4 
addresses the problem of optimal ordering when backorder 
costs are high. This situation corresponds to that of the 


BVP for items which should "never be out of stock". 
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